Low-energy renormalization of the electron dispersion of high-Tc superconductors 
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High-resolution ARPES studies in cuprates have detected low-energy changes in the dispersion 
and absorption of quasi-particles at low temperatures, in particular, in the superconducting state. 
Based on a new 1 /N expansion of the t- J-Holstein model, which includes collective antiferromagnetic 
fluctuations already in leading order, we argue that the observed low-energy structures are mainly 
caused by phonons and not by spin fluctuations, at least, in the optimal and overdoped regime. 
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ARPES studies in cuprates have recently detected changes in the electronic spectral function near the Fermi 
energy at low temperatures and, in particular, in the superconducting state The dispersion of quasi-particles 

changes its slope by about a factor 2 over an energy region of about 50 meV and this change occurs in a rather 
isotropic way around the Fermi surface. These results indicate that the particles near the Fermi surface interact with 
excitations characterized by an energy scale of about 50 meV with a dimensionless coupling constant of about 1. 
Possible candidates for these excitations are phonons or spin fluctuations. In particular, the collective spin excitation 
associated with the resonance peak of neutron scattering has been assumed to cause this new low-energy scale ||] . 

Band structure calculations yield an overall dimensionless coupling constant A of electrons near the Fermi surface 
and phonons between 0.5 and 1.5 [^,0. On the other hand, transport relaxation times indicate a very weak cou- 
pling \tr < 0.3 between electrons and phonons |^ in the optimally doped case. One explanation for this is that 
electronic correlations make the effective electron-phonon coupling strongly momentum and frequency dependent 
suppressing large- momentum scattering . Below we will study whether a similar large reduction of the effec- 
tive electron-phonon coupling by correlations occurs in the electronic self-energy and the renormalization function 
Z which determines the change in the quasi-particle dispersion. Another relevant point is the anisotropy of this 
renormalization. Assuming a bare, momentum-independent electron-phonon coupling the same correlation-induced 
momentum dependence in the effective coupling which suppresses transport scattering rates will make Z in general 
anisotropic. The calculation has to show whether this induced anisotropy is small enough to account for the rather 
isotropic behavior of the slope change observed in experiment. 

Spin fluctuations within the usual 1 /N expansion yield only small contributions to the electronic self-energy and to 
scattering rates at low energies Though the coupling between spin fluctuations and electrons is set by the large 
hopping integral t, the spin fluctuation spectrum at large N is rather structureless and extends in energy over several 
t's. Most spin flips thus are high-energy excitations yielding only small renormalizations at low energies. Below 
we will present results from a different 1/N expansion which favors antiferromagnetic fluctuations and accounts for 
the resonance peak in the spin susceptibility already in leading order. This approach is convenient to make a fair 
comparison of phonon- and spin fluctuations-induced contributions to the electronic self-energy. 

In the following we consider a slightly modified version of the Hamiltonian Eq.(l) of Ref. describing the t-J 
model coupled to dispersionless phonons. First, we assume that the underlying internal symmetry is the symplectic 



group Sp(N/2) and not SU(N). Consequently, we split the internal index p in Eq.(l) of Ref. |11| into an index a 
describing the two spin components and a flavor index fj. running from 1 to N/2. Secondly, we modify the second 
term in Eq.(l) in the following way: 

*J ij, aa' = 1,2 

p,p' = l.-.N fj,,ii' = l...N/2 

In the modified term the sum over the flavor index is carried out independently in each X operator, whereas the sum 
over spins retains the old form. The two terms in Eq.(|l|) are identical in the physical case N=2. At large N they 
differ, however, and the label arrangement in the left and right terms in Eq.(^ favor RVB and antiferromagnetic 
correlations, respectively. Thirdly, we multiply the Hamiltonian with an overall factor 2 so that it coincides with the 
usual Hamiltonian for iV = 2 and t and J have the usual meaning. 

The 11 element of the electronic Green's function matrix in a BCS superconductor has the following form, 

«ti;„ - e(k) - i;(icj„,k) - A2(k)/(iti;„ -f e(k) - S(iti;„,k)) 

u)n = (2n-|-l)7rT is a fermionic Matsubara frequency and T the temperature. A(k) = /S.(cos{kx) —cos{ky)) describes a 
frequency-independent superconducting order parameter with d-wave symmetry which is a rather good approximation 
for the anomalous self-energy at large N [Q. The one-particle energy e(k) is in the leading 0(1) of the 1/N expansion 
given by —2t{cos{kx) + cos(ky)) — At' cos{kx)cos{ky) . t = 1^5/2 and t' = t'f^5/2 are renormalized hopping elements 
with to being bare hopping matrix elements between nearest and second- nearest neighbors on the square lattice, 
respectively. 5 is the doping away from half- filling. In our modified 1/N expansion e(k) does not have any contribution 
from J as long as there is no long-range antiferromagnetic order. E is the diagonal self-energy due to the many-body 
interactions in H. In the leading 0(1/N) S is additive in a purely electronic part Sei and a phononic part, Sp/i, 
renormalized by electronic correlations. We have analyzed all the contributions to Sgj and found that the RPA-like 
term 

^RPA{k) = J(MiHih)G'^"Hk + fci), (3) 
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with the spin vertex 7s(fci) = 1/(1 — ^(ki)a(fci)/2), is the only term which may give rise to interesting structures 
at low energies. In Eq.(^) we combined the Matsubara frequency iujni and the two-dimensional wave vector ki into 
a three-dimensional vector fci. Nc is the number of cells and a(fci) is defined below. The imaginary part of the 
remaining terms of 0(1/N) exhibits a featureless lu'^ law at low energies and acquires structures only on the energy 
scales J or t. This contribution is similar to the total 0(1/N) contribution to the imaginary part of S in the usual 
1 /N expansion depicted in Fig. 5 of Ref. |l^ . 

The phonon contribution Sp/j is given in 0(1/N) by 

Sp.(fc) = ^ E + ^i)? ^^^4— ^7c(fci + k, -fci)7c(fc, fci), (4) 

where ujq and A denote the phonon frequency and the dimensionless electron-phonon coupling constant, respectively, 
and G'-^'' is given by Eq.(^ with S = 0. The phonon vertex 7c reads 

^(kk,)^ ~^-b{ki) + 2a{h)m 

' ^' l + 6(fci)(l + c(fci))-a(fci)d(fci)-|-a(fci)J(ki)/2- ^ ' 

i(k) is equal to — e(k)/((5/2). The four susceptibilites a,b,c,d are the generalizations of the normal state suscepti- 
bilities given in Eqs.(10)-(13) in Ref. (l| to the superconducting state. Omitting the term proportional to J in the 
denominator the expression for 7^ in Eq. (||) is identical with the J ^ limit of the expression obtained in the 1 /N 
expansion used in Ref. j|] . 

Fig. 1 shows the imaginary part of the electronic self-energy in the superconducting state at T = as a function 
of the frequency uj for two momenta Q on the Fermi surface, namely, Q — (1.33, 1.33) — L along the diagonal (left 
panel) and Q — (3.14, 0.26) = X (right panel). The energy unit is an effective nearest neighbor hopping constant t of 
150 meV. t' /t is equal to -0.25 and the value for J, J = 1.11, is chosen such it reproduces the resonance peak in the 
spin susceptibility at the experimental value of about 40 meV. The doping i5 away from half-filling is (S = 0.20. The 
definition of the coupling constant A, Eq.(7) of Ref. is based on the bare coupling and an average density of states 
of l/8t. The dashed and solid lines have been calculated neglecting and including vertex corrections, respectively. 
The dash-dotted lines represent the spin fluctuation contribution, Eq.(H). 




FIG. 1. Imaginary part of the electronic self-energy E(Q,a;) at the points L (left panel) and X (right panel). Solid, dashed, 
and dot-dashed lines correspond to the phonon contributions with and without vertex corrections, and the spin fluctuation 
contribution, respectively. Energies are given in units of t — 150 meV. 



The solid and dashed lines in Fig. 1 show a gap near the Fermi energy bounded by two sharp peaks with an energy 
distance of about 2(A + llIq). The lines fall to zero below uj ^ —2.5 and above u; ~ 6 reflecting the finite electronic 
band-width. Comparing the left and right panel one recognizes that Ep/j is rather isotropic in spite of the fact that 
the vertex corrections depend strongly on momentum and that the superconducting gap is anisotropic. The vertex 
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corrections 7c, which represent one effect of electronic correlations, suppress large momentum transfers in the sum over 
ki in Eq.(^) but enhance small momentum transfers. Comparison of the dashed and solid lines shows that altogether 
the suppression prevails making the solid curves smaller by roughly a factor 2 compared to the dashed ones. 

In order to discuss S_rpa let us first consider the momentum integrated susceptibility at large N, shown 

in the left panel of Fig. 2. In accordance with previous work x"(Q,a;) develops a bound state in a very localized 
region around Q = (tt, tt) with nearly no spectral weight left at higher energies at this momentum. Fig. 2 shows 
that also x"{^) exhibits a sharp structure in the region around 40 meV. It consists of a stronger component at lower 
energy which comes from transitions near (7r,7r) and corresponds to the resonance peak, i.e., an approximate pole 
in the denominator of 7^. The upper weaker peak is mainly due to momenta near the X point and is caused by 
spin flip excitations across the superconducting gap. Experimentally, a two peak structure has been seen in inelastic 
neutron scattering in the momentum integrated spin susceptibility in Y Ba2Cu30Q,5 and Y Ba2Cu30e,'/ (Figs. 5 
and 6 of Ref. Q) but not in x"(Q,cj) for Q = (tt, tt) in rBaaCusOe.T (Fig. 7 of Ref. Q). The higher peak is 
sensitive to temperature, in contrast to the lower peak, and vanishes somewhere above the superconducting transition 
temperature Tc- Guided by the result of our calculation we identify the upper peak in the experimental x"('^) with 
the pseudogap with energies of about 60 and 55 meV in YBa2Cu30e,5 and YBa2Cu30e,i, respectively. Accepting 
this interpretation the pseudogap increases slightly with decreasing doping which is in line with the well-investigated 
behavior of the pseudogap in the Bi-cuprates. Fig. 2 also illustrates that x"(w) extends over a large energy region 
set by the band-width and that most of the spectral weight resides at high energies. Using the exact sum rule one 
finds that the resonance peak exhausts only about 1.5 per cents of the total spectral weight which agrees well with 
experiment |15 1^. 
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FIG. 2. Left panel: Momentum-integrated susceptibility x"(.^ + ^v) using the same parameters as in Fig. 1; right panel: 
renormalization function Z{Q,lj) for interaction with phonons ("Ph") at the point Q = L for T — (solid line), T = ISO-ft' 
(dashed line), and with spin excitations ("RPA") at Q = L and T = (dot-dashed line) and at Q = X and T = (dotted 
line) . 



The dash-dotted lines in Fig. 1 represent the imaginary part of Eijpyi- In contrast to Sp/j SppA is strongly 
anisotropic along the Fermi surface. The low-frequency spin fluctuations which are concentrated around momentum 
transfers of about (tt, tt), can probe the superconducting gap if Q is near the X but not the L point. Therefore 
sharp gap features appear only in the right but not in the left panel. One also notes that the dot-dashed lines are 
substantially smaller in absolute value at low energies and decay slower towards higher energies compared to the case 
of phonon scattering. The reason for this behavior can be inferred from x"(a;), as shown in the left panel of Fig. 2. 
Most of the spectral weight in x"('^) is found at high energies which leads to the slow decay of the dash-dotted curves 
towards high energies. The rather small spectral weight associated with low energies and the resonance peak yields, 
together with J as the coupling constant between electrons and spin fluctuations, only small values in the self-energy 
near the Fermi energy. The neutron resonance peak leads in the calculation of Ref. [|j to much larger effects in S 
and Z. There the employed coupling constant and the frequency-integrated resonance peak are taken as 0.65 eV and 
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13 per cents of the total sum rule, respectively. The corresponding values are in our case J = 0.166eF and 1.5 per 
cents, respectively. Our rather small spectral weight of the resonance peak agrees well with the experiment |l5| , p^ . 
The present 1/N expansion also identifies in a unique way J as the coupling constant between electrons and spin 
excitations. 

The right panel in Fig. 2 shows the renormalization function Z{Q,uj) = 1 — i?e(I](Q,cj) — Y.{Q,0))/uj, where Re 
denotes the real part. The solid line is calculated from the phonon induced self-energy at the point L and at T = 0. 
A very similar curve is obtained for the point X. The solid line deviates substantially from 1 only in a narrow energy 
range around the Fermi energy set by the phonon energy and the superconducting gap assuming a value of 1.76 at 
iv = 0. Defining a renormalized couphng function A(Q) by Z(Q,0) = 1 + A(Q) one finds a rather isotropic A(Q) 
of 0.76. There are two competing effects due to electronic correlations in renormalizing the original bare coupling 
constant A: Vertex corrections suppress T,ph by roughly a factor 2, but density of state effects increase the effective 
A again so that A and A do not differ much. Calculating the related quantitiy Xtr, which determines the coupling 
strength in transport phenomena, one has for the bare values Xtr = A because the electron-phonon interaction was 
assumed to be momentum- independent. The renormalized quantity At,- is, however, equal to 0.24, i.e., it is about 
three times smaller than A. In this case the strong suppression of large- momentum scattering by vertex corrections 
dominates. Such a big difference between A and Xtr is needed to have agreement with experiment: The observed slope 
change in the quasi-particle dispersion by about a factor 2 corresponds to A ~ 1, whereas the transport data suggest 
a three or four times smaller Xtr in the optimally doped case. 

The dashed curve in the right panel of Fig. 2 describes the phonon-induced Z at the L point at T=180K. It 
illustrates that its value at w = is sensitive to temperature and thus would decrease the slope change substantially 
with increasing temperature. The dotted and dash-dotted curves in Fig. 2 represent the spin fluctuation induced 
renormalization factor. It is very anisotropic and in general much smaller than the phonon-induced contribution to Z. 
As shown in Ref. pof the transport quantity Xtr tends generally to be larger than A in the case of spin fluctuations. 
The dotted and dash-dotted curves thus yield an effective Xtr of about 0.2 which is compatible with transport but 
indicates, on the other hand, that the spin fluctuation contribution to Z is not the dominant one. 

The left panel in Fig. 3 shows the evolution of the spectral function at low temperatures in the superconducting state 
as a function of frequency for six momenta approaching the Fermi wave vector in the direction L. The interaction with 
phonons as well as the parameters of Fig. 1 were used in calculating the self-energy. The dashed lines represent the 
spectral function without the self-energy. Going away from the Fermi surface in momentum space the quasi-particle 
peak disperses much weaker than the free peak away from the Fermi energy and looses rapidly spectral weight. The 
weak hump seen near the Fermi surface disperses roughly with the free particle dispersion at somewhat lower energy 
and becomes broad further away from the Fermi surface. Plotting the same as a function of momentum for a fixed 
energy one essentially obtains Lorentzians dispersing like free particles. All these properties are well reflected in 
the experimental data [Q-Q showing that the self-energy is rather isotropic along the Fermi surface but strongly 
frequency-dependent . 
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FIG. 3. Left panel: Spectral function for phonon scattering for 6 momenta approaching the Fermi vector along the 
(tt, tt)— direction, dashed curves are obtained neglecting self-energy effects; right panel: imaginary part of n(Q,a;) for three 
temperatures, using parameter values as in Fig. 1. 



From this one concludes again that the spin fluctuation contribution to Z cannot dominate in agreement with the 
right panel in Fig. 2. 

The right panel in Fig. 3 shows the imaginary part of the total self-energy 11 associated with G, i.e., n(zw„,k) — 
—G~^{iujn, k)-|-iLj„ — e(k). The curves in this Figure can directly be compared with the corresponding curves extracted 
from ARPES data. Fig. 2a of Ref. |17|. In the calculation we used for all temperatures the same T=0 gap A which 
approximately is found in the analysis of the data jl^]. There is good agreement on a qualitative level. Quantitatively, 
there are, however, several discrepancies: in order to reproduce the absolute value of the hump (our energy unit t 
corresponds to 150 meV) A should be increased by about a factor 2 which would tend to yield too large slope changes 
in the quasiparticle dispersion at the point L. At low temperatures the gap between the hump and zero energy is in our 
case much more pronounced than in the experiment. One reason for this discrepancy is our use of just one phonon 
frequency instead of the true, rather broad phonon density of states. It is also possible that the spin fluctuation 
contribution cannot be totally neglected at the point X and enhance the dump-dip feature. 

The situation is less clear in the underdoped region. Experimentally, the momentum-integrated susceptibilities, if 
integrated between and 65 meV, are typically one order of magnitude or more larger than in the optimally doped 
case. For instance, YBa2Cu30e,r exhausts in this region about 17 (10) per cents of the sum rule at T=12K (200K) 
10 whereas the resonance peak does only 1 per cent ||l^,0 . Using the experimental susceptibility in Eq. (^) instead 
of the calculated large N susceptibility we find in the above case for Q = L Z{Q,0) = 1.47(1.27) and for Q = X 
^(QjO) — 3.30(2.31) at T=12K (200K). Because Xtr > A in the case of spin fluctuations the resulting values for Xtr 
seem to be too large to be compatible with the resistivity curves [|l8|. Also the large anisotropy in the slope change 
of the quasi-particle dispersion, implied by these values, seems not to be supported by the available experiments. 

Using a 1/N expansion which takes into account the magnetic resonance peak already in leading order we have 
calculated the contributions of 0(1/N) to the electronic self-energy for a t-J-Holstein model. We find that the recently 
observed low-energy renormalization of quasi-particles in the superconducting state of cuprates can well accounted 
for taking a phonon energy of 35 meV and a dimensionless coupling constant of about 0.75 — 1. In calculating the 
spin fiuctuation part to the self-energy our approach does not contain any adjustable parameter. Due to the rather 
small value of J, which acts as the coupling constant, and due to the structureless x" large N we find that this 
contribution is rather small and strongly anisotropic along the Fermi surface. This means that it does not play the 
dominant role in the observed low-frequency renormalization of quasi-particles, at least, in the optimal and overdoped 
regime. 
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